1 OPEN ARCHITECTURE GRADIENT COIL SET FOR MAGNETIC 

2 RESONANCE IMAGING APPARATUS 
3 

4 

5 CROSS-REFERENCE TO RELATED APPLICATION 
6 

7 This application claims the benefit of U.S. provisional patent apphcation Serial No. 

8 60/270,960 filed Febraary 22, 2001. 
9 

10 BACKGROUND OF THE INVENTION 
11 

12 The present invention relates to magnetic resonance imaging and, in particular, to 

13 gradient coils used therewith. 

, 14 A Magnetic Resonance Scanner where the direction of the main static field is orthogonal 

^15 to the pole surfkce and to the magnet's main symmetry plane is commonly called an Open 

i;|6 Magnet System. 

;; i 1 Previous methods for production of linear magnetic gradients in such MRI systems 



i,a 8 consist of winding discrete coils, in a bunched or distributed fashion on an electrically insulating 

ri9 hollow right cylinder coil-form, and driving said coil with a limited voltage. Conventional 

llilo bunched coil designs are the Maxwell and Modified Maxwell Pair for the Z-gradient coil 

(parallel to the main magnetic filed), and the Golay or Modified Golay (multi-arc) Saddle coils 

i^i|2 for the X/Y gradient coils. Previous methods consisted of iteratively placing elemental loops or 

' %3 arcs on the cylindrical coil form until the desired gradient strength, gradient uniformity and 

24 inductance were achieved. In general, these types of designs constituted as method referred as a 

25 "forward approach" whereby a set of initial coil positions (i.e. initial current distribution) were 

26 defined. The gradient field, the energy/inductance were calculated and if they did not fall within 

27 specified design criteria, the positions of the coils will be shifted (statistically or otherwise) and 

28 the new outcome will be re-evaluated. This iterative procedure will continue until a suitable 

29 solution that satisfies all the design criteria is found. 

30 An altemative approach in the design of gradient coils is referred as an "inverse 

3 1 approach". According to this methodology, a set of predetermined constraint points for the 

32 gradient magnetic field is chosen inside the desired imaging volume. Via a Lagrange 

33 minimization technique involving the stored magnetic energy or power of the structure, the 
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1 appropriate continuous current distribution that generates the predetermined gradient field inside 

2 the imaging volume is derived. Using a stream function discretization technique, the continuous 

3 current distribution of the coil can be approximated by a set of discrete current loops that are 

4 connected in series and share the same constant current. In order to ensure that the discretization 

5 process is valid, the gradient field is regenerated using the Biot-Savart law to the discrete current 

6 distribution and compare the results with the ones of the analytic approach. If the error is 

7 relatively small, the discrete current distribution is accepted as a valid one. 

8 Open magnet geometries with vertically directed fields have been a successful presence 

9 in the MRI arena due to the inherent open architecture scheme that promotes patient comfort and 

10 diminishes patient claustrophobia. Recently, vertical field systems with higher main field 

1 1 strengths (0.5 T to 1 .0 T) were also introduced in an attempt to improve image quality, as well as 
^12 to expand the spectrum of clinical applications, especially fast imaging techniques not available 

on lower field strength MR scanners. The necessity for utilizing such imaging techniques is 

•^.|4 strongly dependent upon the availability of a suitable gradient system capable of delivering high 

li'fs gradient field strengths with a significant improvement to the gradient coiFs slew rate. For open 

!,|6 magnet geometries with vertically directed fields, biplanar gradient coils have been the geometry 

"' t? of choice, since there are conformed to the main magnet's geometric shape and assist in 

:;l8 improving patient comfort and reducing patient claustrophobia. There is however a drawback in 

1 1 9 regards to the biplanar designs. Because of the demand that the coil must have enough gap to 

;,;|0 allow for clear access to the human torso, the distance between the two planes is excessively 

i'^ll high. Increasing the gap between the two planes reduces the efficiency of the biplanar gradient in 

22 a level. Thus for a whole body gradient coil set, it is very difficult to generate gradient fields of 

23 the order of 50-60 mT/m with slew rates approaching 1000 T/m/sec. An alternative solution to 

24 this issue is the introduction of xmiplanar gradient coil designs, which offer superb gradient field 

25 characteristics over a very limited imaging volume. Uniplanar gradients are mainly utilized on 

26 local imaging applications and are not recommended for apphcations with sizable FOV such as 

27 the one covering the human head. One viable alternative for performing ultra fast imaging of the 

28 human head for a vertical field system is the use of insertable biplanar gradients in a very close 

29 proximity to the human head. Due to the reduction of the gap between the planes the efficiency 

30 of the coil is significantly improved, and such a design is capable of generating high peak and 

31 high uniformity gradient fields over a volume similar to the size of the human head (25 cm 



2 



1 DSV). But such a design is prohibitive because it violates the two principles of which the open 

2 magnet designs are based on. Placing the two planes over tiie human head and especially over 

3 the patient's eyes, it compromises the patient's comfort and enhances patient anxiety. 

4 Gh-adient coil sets in common use make serious tradeoffs between coil performance and 

5 maintaining the "open," non-claustrophobic configuration, desired in open MRI systems. 
6 

7 SUMMARY OF THE INVENTION 
8 

9 A MRI gradient coil set includes a uniplanarZ-gradient coil, a biplanarX-gradient coil, 

10 and a biplanar Y-gradient coil. The coil set provides an open Z-axis face. 

11 

12 BRIEF DESCRIPTION OF THE DRAWINGS 

s 

FIG. 1 is a perspective view of a possible gradient coil set according to the invention. 

P FIG. 2a is a table showing on-axis linearity. 

I FIG. 2b is a table showing off-axis uniformity. 

% FIG. 3 is a table showing on-axis field behavior. 

i;I|8 FIG. 4 is a table showing on-axis field behavior. 

;ip FIG. 5 is a gradient pattern. 

ilO FIG, 6 is a gradient pattern. 

f|l FIG. 7a is a gradient pattern. 

22 FIG. 7b is a gradient pattern. 

23 FIG. 8 is a table showing on-axis field behavior. 

24 FIG. 9 is a gradient pattern. 

25 FIG. 1 0 is a gradient pattern. 

26 FIG. 1 1 is a gradient pattern. 

27 FIG. 12 is a gradient pattern. 

28 FIG. 1 3 is a gradient pattern. 

29 FIG. 14 is a table showing on-axis behavior. 

30 FIG. 15 is a gradient pattern. 

3 1 FIG. 1 6a is a gradient pattern. 
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1 FIG. 16b is a gradient pattern. 

2 FIG. 17a is a gradient pattern. 

3 FIG. 17b is a gradient pattern. 

4 FIG. 1 8 is a perspective view of a possible gradient coil set according to the invention 

5 with shielded biplanar coils. 

6 FIG. 19 is a perspective view of a possible gradient coil set according to the invention 

7 with unshielded biplanar coils. 

8 FIG. 20 is a perspective view of a possible gradient coil set according to the invention 

9 with shielded biplanar coils and incorporating radio frequency coils. 

10 FIG. 21 is a perspective view of a possible gradient coil set according to the invention 

1 1 with unshielded biplanar coils and incorporating radio frequency coils. 

12 FIG. 22 is a perspective view of a possible gradient coil set according to the invention in 
i;||3 an array configuration. 

!'jC4 FIG. 23 is a perspective view of a possible gradient coil set according to the invention 

Jjis incorporating knee and upper thigh frequency coils. 

i FIG. 24 is a perspective view of a possible gradient coil set according to the invention 

T7 incorporating knee and foot radio frequency coils. 

m 

i|9 DESCRIPTION OF THE PREFFERRED EMBODIMENTS 




fil In this invention, a new alternative design for an insertable open architecture gradient coil 

22 set is given which is focused on enhancing patient comfort and abolishing patient anxiety, while 

23 delivering a gradient field with high peak strength and improved slew rate. According to this 

24 architecture, the Z gradient coil is a uniplanar design that is placed in parallel with the main 

25 magnet's poles. In order to form a complete gradient set, the two transverse QC, Y) gradients are 

26 generated by a pair of biplanar designs which are placed perpendicular to the magnets poles, but 

27 they may be physically mounted on the same side with the uniplanar Z gradient. The gradient 

28 field generated from the two transverse coils exhibits similar peak and slew rate performances as 

29 the uniplanar Z gradient coil over a comparable imaging volume. When this coil is designated for 

30 head imaging applications, p^abolic cut-offs on the two side transverse biplanar coils may be 

31 implemented in a symmetric fashion. Such an insertable coil design does not disrupt the vision of 



4 



1 the human subject and thus makes it more patient friendly without imposing restriction on the 

2 field strength and quaUty. The proposed design can also be implemented in a self-shielded 

3 configuration for reduction of induced eddy current effects on the vicinity of the main magnet's 

4 poles. In addition, the proposed design can be paired with an additional set of gradient coils 

5 (shielded or not) in an array configuration with a capability of extending the imaging coverage of 

6 the gradient coil system, while improving patient comfort and diminishing patient 

7 claustrophobia. 
8 

9 Theory 
10 

1 1 FIG* 1, shows the geometrical shape for a self-shielded open architecture gradient set. 

12 The gap between the two primary planes on the biplanar section of the coil is denoted as 2a while 
If the distance between the shielding planes of the biplanar section of the coil is denoted as 2b. For 

the uniplanar section of the coil, the distance between the plane section of the coil and the 
geometric center of the coil is denoted as a', while the distance between the geometric center of 

])S the coil and shielding part of the biplanar section as b'. The schematic representation of a head 

'1*7 Transmit/Receive Radiofrequency coil is also included. 

\ IB 

119 Axial (Z) Uniplanar Gradient 

I 

Wl Let us assume that the position of the plane for the primary gradient coil is at z=-a' from 

22 the geometric center of the coil, while the position of the shielding coil is z=-b'from the center of 

23 the coil. The current density distribution for this design is confined on the xy plane and has the 

24 expression: 

J ix, y) = [J^ (x, y)x + Jy (x, y)y] (1) 

25 
26 
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1 In addition the expression for the vector potential is given by: 
2 



1 



>-<i (2) 



3 
4 



And for the region outside the planes z<-b', the expression for the vector potential is: 



4 =A J 



1 



where J(x,y) (a,P) represent the double Fourier Transform of J(x,y) (x,y), respectively. Since there 
is no charge free in the vicinity of the gradient coils, the application of the current continuity 
equation (V»J=0) leads to a relationship between the components of the current density as: 



a 



i:;38 JJa,p) = —-J,(a,P) 



(6) 



From equations (2)-(5) and equation (6), the expression for the vector potential at the location of 
the two planes is given by: 



1 



1 



forz = -a' 



(7) 



a 

7. 



forz = -a' (8) 



'j/ia,P)+J,-^{a,P)^forz=-b' 



4 =^ j £dc^^-jJ==e^--^l-^J^^ 



(9) 
(10) 
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12 The expression for the magnetic field is: 



6 



gs 
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11 

12 



[l-e-'^<*'-'''^]/orz>-a' (11) 



y5) 



1 
2 

3 and since the magnetic field is symmetric along the x and y directions and monotonically 

4 increased along the z direction, equation (1 1) becomes: 



with 



Mo 



J £ia^A^^^^cos^)cos(^K^^^"V(«.y^) 



1-e 



(12) 



j;''\a,J^=-e^^^~'''^j;''\a,fi)ishielding conditionB^ =Oforz>-b') (13) 
and 

'^/icc,^=i^x"'i(x,jS)^j/(a,/3)=4j dxdy cos^) sm(j9^)j;''' ix,y) (14) 



Fuithermore, the expression of the stored magnetic energy is: 



(15) 



In order to find the optimum continuous current density, which will generate a gradient magnetic 
field in predetermined locations inside the imaging volume, can be extracted via a Lagrange 
optimization technique, According to this technique, a functional expression, which is a 
quadratic function of the current density is constructed. An example of such an expression 
includes which includes the stored magnetic energy and the magnetic field, and can have the 
form: 



(16) 



13 where are the Lagrange multipliers and Bzsc(rj) are the constrained values of the magnetic field 

14 inside the desired imaging volume. Minimizing F with respect to Jx'^ , a liaear matrix equation 

1 5 for Jv'^ is obtained: 
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j/ = -fij^^j cosiccXj)cos(ffyj)e-^''^''''' (17) 
1 and the Lagrange multipliers are determined via the constrained equation 



+C0 



2 

3 For the case of a non-shielded nniplanar axial gradient set, the above expressions are the 

4 same at the limit v^^hen b'->oo. 
5 

4 From the continuous current distribution and using employing the stream function technique, a 

; j7 discrete set of loops that closely approximates the continuous current density pattern is 

3 generated. In order to verify the accuracy of the discretization, the gradient field of the coil is re- 

4 

\ ,3 evaluated by applying the Biot-Savart law to the discrete current distribution. 



m 

iii|2 Transverse X or Y Uniplanar Gradient 

lis 

i;:3 

f|4 The expression for the magnetic field is fix)m equation (1 1): 

=-dA f rdodfi^[^^'-^^e-^^'*''j,-'(a,^ 

\\-e'-^'''-'''^forz>-a' 

15 
16 

17 and since the magnetic field is symmetric along the x and z directions and antisymmetric along 

18 the y direction, equation (11) becomes: 
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(12«) 



J/*' (a, y5) = j^-"' (a, fi){shielding condition : B, = z > -b') (13fl) 

J/"' {a,j3) = 4j ^dxdy sin(ca) sinCy^j) J/"' (jc, 3;) (14a) 



1 
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Furthermore, the expression of the stored magnetic energy is: 



'.=Ud^x A.J=^f\dadp^^^^\j;\aM 



(15a) 



(3 
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In order to find the optimum continuous current density, which will generate a gradient magnetic 
field in predetermined locations inside the imaging volume, can be extracted via a Lagrange 
optimization technique, According to this technique, a fimctional expression, which is a 
quadratic fimction of the current density is constructed. An example of such an expression 
includes which includes the stored magnetic energy and the magnetic field, and can have the 
form: 



(16) 



where are the Lagrange multiphers and B2sc(rj) are the constrained values of the magnetic field 
inside the desired imaging volume. Minimizing F with respect to Jx"* , a linear matrix equation 
for Jx"* is obtained: 



j/ = -^Y, sin(cix^ ) cos(y^^ w) ^^^^^ 



12 and the Lagrange multiphers are determined via the constrained equation : 
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tc^^j =B.sc, =^C, = ^ J ]dadfismiax,)cosi^,)e-^'''^'''' 
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1 For the case of a non-shielded uniplanar axial gradient set, the above expressions are the 

2 same at the limit when b'->oo. 

3 

4 From the continuous current distribution and using employing the stream function technique, a 

5 discrete set of loops that closely approximates the continuous current density pattern is 

6 generated. In order to verify the accuracy of the discretization, the gradient field of the coil is re- 

7 evaluated by applying the Biot-Savart law to the discrete current distribution. 
8 

9 Phased Array considerations: Mutual energy between two coils. 
10 

1 1 Let us assume a second planar set with planes locations of z— a ' and z=-b\ If the second 

1,12 set is shifted along the x or y directions or both by an amount of xo and yo, respectively, the 



!;6 analytical expression of the current density distribution is given by: 

% 

sil? Considering for simplicity of the evaluations, that the second set is shifted only along the y 

f ij8 direction by an amount of yo. In this situation the mutual energy between these two modules is 

I II 9 given by: 

m 

21 -i P 



j_^-2V^(6'-«.'-«0)j (20) 



22 

23 If both planes represent the same gradient coil axis, equation (20) becomes 



+w / 2 I 

W.^=-^[ \da dfi ^ / cos.{fy,)j;''\a,p)Jr'\a,fi) 
24 4;r _i fi 



(21) 
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The optimum phased array configuration occurs at the yo points where the value of the mutual 
inductance between the two sets is minimized. 



Transverse (X, Y) Biplanar Gradient Coils 



Let us assume that for the self-shielded biplanar design, the gap between the two planes of the 
primary gradient coil is y=2a, while the gap between the planes of the shielding coil is z=2b. 
The cuirent density distribution for either the primary gradient coil or the secondary gradient coil 

[ j;' (x, z)x + (jc, z)z] 5{y-b)+ {j;' (x, z)x + j;' (x, z)z] 6{y + b) (22) 
is confined on the xz plane. The analytical expression for the current density of the coil is: 

In addition the expression for the vector potential is given by: 



where J(x,z) (a^p) represent the double Fourier Transform of J(x,z)(x,z), respectively. Since there 



for\y^a (23) 



for\y^b 



is no charge free in the vicinity of the gradient coils, the appUcation of the current continuity 
equation (V«J=0) leads to a relationship between the components of the current daisity as: 



a 



(25) 



From equations (23)-(24) and equation (25), the expression for the vector potaitial at the 
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1 location of the two planes is given by: 



/;(a,;9)+e-^^^Vr(«,;6)+ 



fory=a 



{26) 



( a' 



fory=a (27) 



/ory=-fl (2^ 



^2 I P) 



fory=-a (29) 



/or>'=6 (3(^ 



/or>'=& (31) 



^3 
4 



Mo 



lax+ipz 



e^<*^''V/(a,)8)+e-^'^<'*">J/''(a,y?) + 



fory = -b (32) 
e^<*^">//(a,^) + e-''^<*^'">y,-"(a,A) +' 



fory^-b (33) 



The expression for the magnetic field is: 

B-^(*-^>y/(a,;ff)-e-^^*^^>j/(«,;0) 



B. 



e-^^-*V/(a,y?)+e^<*"^V/\a,;ff) 



/orlylSa (34) 



/or|y^6 (35) 



5 
6 
7 



Transverse Y Gradient Coil 
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1 For the Y gradient coil, the gradient magnetic field is symmetric along the x and z directions 

2 while it is asymmetric along the y direction. In this case the relationship between the current 

3 densities on the adjacent planes is: 

j:{a,P) = j;''{a,P) and J ^\a,p) = J ;\a,P) 4 

5 

6 The expression of the magnetic field firom equations (34) and (35) becomes: 



/;(a,;^(e^'^^^ -e^<^^>) 



Jcoc+i/k 



//(a,/7)(-^-^<^> 



for\y\<a (36) 
for\y^b 



(37) 



. -7 



''"3 



m 



The gradient field shielding condition in the region outside both planes demands that Bz=0 for 
|y|>b. Equation (37) leads to: 



j;"' {a,f3) = -e-^'^'''-'''^j;"\a,P) (38) 
whichleadsto : 



= --^ J £^dcdficos(coc) cosme-^" sinh(yy[aU^)j; (a, J3^ -g-^^^"-") (39) 



and 



11 



//(a,>ff) = 4j ^dxdz cos(ccc)cos(j^)J°(x,z) (40) 
Employing the current density's symmetry conditions, the shielding condition and the continuity 



12 
13 
14 
15 
16 
17 



equation, the expression for the stored magnetic energy is given by: 



l-e 



(41) 



In order to find the optimum continuous current density, which will generate a gradient magnetic 
field in predetermined locations inside the imaging volume, can be extracted via a Lagrange 
optimization technique. According to this technique, a functional expression, which is a 
quadratic function of the current density is constructed. An example of such an expression 
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1 includes which includes the stored magnetic energy and the magnetic field, and can have the 

2 form: 



F(J/)^W^ -•f^ZjlB,(rj)~B^Arj)^ (42) 

3 where Xj are the Lagrange multipliers and B^scCrj) are the constrained values of the magnetic field 

4 inside the desired imaging volume. Minimizing F with respect to Jx^ , a hnear matrix equation for 

5 Jx^ is obtained: 

// = __1L=— ll^^L^f cos(aXj)cosi/^j)sivh(yj^a' (43) 

6 and the Lagrange multipliers are determined via the constrained equation : 



{dadflcosi ax,)cos( fiz,)smK y,^a' + fi') ^ ~ ' , (44 ) 

71 ' ^ ^ ^a^ + fi^ cosh( a^a^ ^ P^) 

sinh(( b - a)'yJa^~V~p^)Y^ ;i ^ cos( a;Cj)cos( fiz j)smh{ y j-Ja^ fi^ ) 

i| 

^% And the expression for the stored energy is: 



w. = 



/?2 -b^a'*p' 

= £^f [dad/3 -i==£= — , sinh(( b-a)^Ja' + fi') (45) 

j^f; cos( ax^)cos( j5z,)sinh( >>^V«' + 

11 

12 

13 For the case of a non-shielded biplanar transverse gradient set, the above expressions are 

14 the same at the limit ivhen b->oo. 
15 
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1 From the continuous current distribution and using employing the stream function technique, a 

2 discrete set of loops that closely approximates the continuous current density pattern is 
generated. In order to verify the accuracy of the discretization, the gradient field of the coil is re- 
evaluated by applying the Biot-Savart law to the discrete current distribution. 



6 Transverse X Gradient Coil 

7 

8 For the X gradient coil, the gradient magnetic field is symmetric along the y and z directions 

9 while it is asymmetric along the x direction. In this case the relationship between the current 

1 0 densities on the adjacent planes is: 

11 J/(a,/?) = -J/^(a,y5) and J^\a,fi)^-j;\a,P) 

12 The expression of the magnetic field from equations (34) and (35) becomes: 



m 

"4 



'■1a 

r|4 
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17 



B=- 



for\y\<a 



(46) 



B=- 



8;r 



/;(a,y&)(-^-^''2cosh^V«' + 1^' ))- 



for\y)^b (47) 



The gradient field shielding condition in the region outside both planes demands that Bz=0 for 
|y|>b. Equation (47) leads to: 



J:\a,/J) =-e-^^*'-"V/'''(«,y0) 
whichleadsto: 



(48) 



B,=-^l £ia;y3sin(a^)cosC&)e-^"cosh6;Va' +>5')>>/(a,y^l-e"'^ 



and 



j;{a,P)=ij;{a,P)^j;{a,^=A\ |°<£c</z sin(ca)cos(^)//(x,z) (50) 
Employing the cmrent density's symmetry conditions, the shielding condition and the continuity 



W.^\\d^x A.J=^^f[dadp^i^\J:iaM[l~e-^^^ 
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1 equation, the expression for the stored magnetic energy is given by: 
2 

3 In order to find the optimum continuous current density, which will generate a gradient magnetic 

4 field in predetermined locations inside the imaging volume, can be extracted via a Lagrange 

5 optimization technique. According to this technique, a functional expression, which is a 

6 quadratic function of the current density is constructed. An example of such an expression 

7 includes which includes the stored magnetic energy and the magnetic field, and can have the 

8 form: 



F{J/)-W„ -YXj[BArj)-B,,,(rj)] (42) 



9 where Xj are the Lagrange multipliers and Bzsc(rj) are the constrained values of the magnetic field 
Ip inside the desired imaging volume. Minimizing F with respect to Jx^ , a linear matrix equation for 
O Jv^ is obtained: 



m 
u 



j; = ^1—^ "/^ Y^j sin(ar,)cos(A,)cosh(7, + ) (52) 
and the Lagrange multipliers are determined via the constrained equation : 



m 



+00 t /? 

r -t^i r^/arf>ffsin(aXj)cos(>&^.)cosh(yy-Ja^ 
14 

15 For the case of a non-shielded biplanar transverse gradient set, the above expressions are 

16 the same at the limit when b->oo. 

17 

1 8 From the continuous current distribution and using employing the stream function technique, a 

1 9 discrete set of loops that closely approximates the continuous current density pattern is 
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generated. In order to verify the accuracy of the discretization, the gradient field of the coil is re- 
evaluated by applying the Biot-Savart law to the discrete current distribution. 



Transverse Z Gradient Coil 



For the Z gradient coil, the gradient magnetic field is symmetric along the y and x directions 
while it is asymmetric along the z direction. In this case the relationship between the current 
densities on the adjacent planes is: 

j:{a,P) = -J:\a,^) and J,\a,P) = -j;\a,P) (54) 

The expression of the magnetic field fi-om equations (34) and (35) becomes: 



5 =- 



.A. 

8^^ 



+ 



for\y^a (55) 



j;(«,y^(-e^"2coshe'V«' +/))- 
j/(a,y^(-ne*^2cosh^V«' +/)) 



for\y\>b (56) 



The gradient field shielding condition in the region outside both planes demands that Bz=0 for 
|y|>b. Equation (56) leads to: 



y/(a,;^=^-^<*'-''V/(a,)8) (57) 
whichleadsto : 



B = 



J ^dodfizQ^) smi^y^^" cosh(vVa' + y^' )^/(a,yS) i-g'^^/^^*-') 



(58) 



and 



j;{a,fi)^ij;{a,j3)^j;{a,p)^^i\ [dxdz cos(at)sinOa) J/(jc,z) (59) 
Employing the current density's symmetry conditions, the shielding condition and the continuity 

equation, the expression for the stored magnetic energy is given by: 
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1 In order to find the optimum continuous current density, which will generate a gradient magnetic 

2 field in predetermined locations inside the imaging volume, can be extracted via a Lagrange 



Q 



-^^jlBArj)-B,,Arj)] (42) 

3 optimization technique, According to this technique, a functional expression, which is a 

4 quadratic function of the current density is constructed. An example of such an expression 

5 includes which includes the stored magnetic energy and the magnetic field, and can have the 

6 form: 

7 where 7^ are the Lagrange multiphers and BzscCij) are the constrained values of the magnetic field 

8 inside the desired imaging volume. Minimizing F with respect to Jx^ , a linear matrix equation for 

9 Jv^ is obtained: 



and the Lagrange multipliers are determined via the constrained equation : 



Q == — J Jdadficos{aXi)sm(^•)coshiyi^Ja^ + fi^ ^ 

(1-^-^^^) § 



Q 1 ^ iyi. cos(6a,)sin(Ay)cosh(j; Ja' 
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12 

13 For the case of a non-shielded biplanar transverse gradient set, the above expressions are 

14 the same at the limit when b-^oo. 

15 

16 From the continuous current distribution and using employing the stream function technique, a 

17 discrete set of loops that closely approximates the continuous current density pattern is 

1 8 generated, hi order to verify the accuracy of the discretization, the gradient field of the coil is re- 

19 evaluated by applying the Biot-Savart law to the discrete current distribution. 
20 
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Mutual Inductance evaluations for Biplanar Open Architecture coils 

Let us assume a second biplanar gradient set with plane gaps for the primary and secondary coils 

at locations of z=2a and z=2b , respectively . If the second set is shifted along the x or z 
directions or both by an amount of xo and zq, respectively, the analytical expression of the current 
density distribution is given by: 



Considering for simplicity of the evaluations, that the second set is shifted only along the z 
direction by an amount of zo. In this situation the mutual energy between these two modules is 
given by: 



, = - \d'x . r^'*' =^ f Ida dfi ^^'t^'e '^^' 

2 J An' J i ^ 



+ 



(63) 



Transverse Y gradient Coil 

Considering the current density symmetries amongst the planes of the self-shielded biplanar 
gradient coil design, as well as that the mutual energy is evaluated over two gradient coil 
representing the same axis, equation (63) has the form: 
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(64) 
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The optimum phased array configuration occurs at the Zo points where the value of the mutual 
inductance between the two sets is minimized. 

X and Z gradient Coils 

Considering the current density symmetries amongst the planes of the self-shielded biplanar 
gradient coil design, as well as that the mutual energy is evaluated over two gradient coil 
representing the same axis, equation (63) has the form: 



+00 / 2 , o2 

^0 J \dadp ^ J cos{fiz,)j:ia,P)J:\a,p) 



An 



(65) 



The optimum phased array configuration occurs at the zo points where the value of the mutual 
inductance between the two sets is minimized. 

Design and Results 

l)Non-Shielded Biplanar Gradient Coils 

i) Transverse Y biplanar gradient coil: 

For the non-shielded transverse biplanar gradient coils, where the gap between the two 
planes is chosen to be 2a=288mm, the constraint points for the minimization algoritiim are 
shown in Table 1 . The gradient field at the center of the structure is demanded to be 40 mT/m 
with 10% on axis linearity at a distance of + 130mm ]from the center of the gradient and 20% of f 
axis uniformity inside a 25 cm Diameter Spherical Volume (DSV). 
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1 



Constraints for a Non-Shielded Biplanar Y gradient coU 


Xj in meters 


Yj in meters 


Zj in meters 


Bzscj in Tesla 


0.000 


0.001 


0.000 


0.0000400 


0.000 


0.130 


0.000 


0.0046800 


0.125 


0.001 


0.000 


0.0000320 


0.000 


0.001 


0.125 


0.0000320 



2 Table 1 

3 

4 The application of the Lagrange minimization technique to the set of constraints of the 

5 Table 1, generates a continuous current distribution for primary coil of the Transverse Non- 

^ 6 Shielded Y biplanar coil. With the assist of the stream function technique, the continuous current 

i;I densities for the two coils can be approximated by a 10 discrete loops with a common current of 

iiji 287.6 Amps. Table 2 illustrates all the vital characteristics needed for the engineering and 

^ manufacturing phase of such a gradient coil. In addition figure 2 shows the on axis field behavior 

l# ( on-axis linearity is represented by Figure 2a, off-axis uniformity is represented by Figure 2b )of 

n the coil as it has been evaluated by applying the Biot-Savart law to the coil's discrete current 

11 pattern. 

m 



Non-Shielded Y Biplanar Gradient Coil 


Property 




Gap 


288 mm 


Gradient Strength at mT/m/lOOA 


13.9mT/m/100A 


Gradient Strength at 500 A 


69.5 mT/m 


450 A 


62.55 mT/m 


330 A 


45.87 mT/m 


Number of Discrete Loops 


10 


Cu thick, between adjacent loops 


> 12.5mm 


Total hiductance (cable included) 


95|jH 


Total Resistance(Cable included) 


SStoQ 
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Rise Time at 2000V/500A 


24 usee 


700V/450A 


70 \isQC 




82 usee 


Slew Rate at 2000V/500A 


2865 T/m/sec 


700V/450 A 


880 T/m/sec 


400V/330A 


550 T/m/sec 


Duty Cycle (RMS Gradient) 




(500 A ref. point) at 4500 kW 


72% of the peak gradient(49.68mT/m) 


at 2000 kW 


49% of the peak gradient(33.81mT/m) 


at 1000 kW 


34% of the peak gradient(23.46mT/m) 



1 Table 2 

Q il) Transverse X biplanar gradient coil: 

IB? 

t$ For the non-shielded transverse biplanar gradient coils, where the gap between the two 

■"I 

planes is chosen to be 2a=284mm, the constraint points for the minimization algorithm are 
i shown in Table 3, The gradient field at the center of the structure is demanded to be 40 mT/m 

7 with 1 0% on axis linearity at a distance of ± 130mm from the center of the gradient and 20% of f 

•''1 

axis uniformity inside a 25 cm Diameter Spherical Volume (DSV). 



Constraints for a Non-Shielded Biplanar X gradient coil 


Xj in meters 


Yj in meters 


Zj in meters 


Bzscj in Tesla 


0.001 


0.000 


0.000 


0.0000400 


0.130 


0.000 


0.000 


0.0046800 


0.001 


0.125 


0.000 


0.0000320 


0.001 


0.000 


0.125 


0.0000320 



10 Tables 
11 

12 The application of the Lagrange minimization technique to the set of constraints of the 

13 Table 3, generates a continuous current distribution for primary coil of the Transverse Non- 
14 Shielded X biplanar coil. With the assist of the stream function technique, the continuous current 
15 densities for the two coils can be approximated by a 10 discrete loops with a common current of 

22 



282.9 Amps). Table 4 illustrates all the vital characteristics needed for the engineering and 
manufacturing phase of such a gradient coil. In addition figure 3 shows the on axis field behavior 
of the coil as it has been evaluated by applying the Biot-Savart law to the coil's discrete current 
pattern. 



Non-Shielded X Biplanar Gradient Coil 


Property 




Gap 


284 mm 


Gradient Strength at mT/m/lOOA 


15.62 mT/m/lOOA 


Gradient Strength at 500 A 
450 A 

330 A 


78.1 mT/m 

70.2 mT/m 
51.5 mT/m 


Number of Discrete Loops 


10 


Cu thick, between adjacent loops 


> 12.5mm 


Total Inductance (cable included) 


177hH 


Total Resistance(Cable included) 


59 mQ 


Rise Time at 2000V/500A 
700V/450A 


45 y.sec 
118 \isec 
154 ^sec 


Linear Slew Rate at 2000V/500A 
700V/450A 
400V/330A 


1738 T/m/sec 
593 T/m/sec 
335 T/m/sec 


Duty Cycle (RMS Gradient) 
(500 A ref. point) at 4500 kW 
at 2000 kW 
at 1000 kW 


68% of the peak gradient(53.11mT/m) 
45% of the peak gradient(35.37mT/m) 
32% of the peak gradient(25.0mT/m) 



Table 4 



iii) Axial Z biplanar gradient coil: 

For the non-shielded Axial biplanar gradient coils, where the gap between the two planes 
is chosen to be 2a=292nim, the constraint points for the minimization algorithm are shown in 
Table 5. The gradient field at the center of the structure is demanded to be 40 mT/m with 10% on 
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1 axis linearity at a distance of ± BOnun from the center of the gradient and 20% of f axis 

2 uniformity inside a 25 cm Diameter Spherical Volume (DSV). 
3 



Constraints for a Non-Shielded Biplanar Z gradient coil 


Xj in meters 


Yj in meters 


Zj in meters 


Bzscj in Tesia 


0.000 


0.000 


0.001 


0.0000400 


0.000 


0.000 


0.130 


0.0046800 


0.125 


0.000 


0.001 


0.0000325 


0.000 


0.125 


0.001 


0.0000325 



4 Table 5 
5 

jS The application of the Lagrange minimization technique to the set of constraints of the 

ffi Table 5, generates a continuous current distribution for primary coil of the Axial Non-Shielded Z 

ill biplanar coil. With the assist of the stream function technique, the continuous current densities 
for the two coils can be approximated by a 10 discrete loops with a common current of 298.15 

W Amps. Table 6 illustrates all the vital characteristics needed for the engineering and 

I i manufacturing phase of such a gradient coil. In addition figure 4 shows the on axis field behavior 
fi of the coil as it has been evaluated by applying the Biot-Savart law to the coil's discrete current 

II pattern. 



Non-Shielded Z Biplanar Gradient Coil 


Property 




Gap 


292 mm 


Gradient Strength at mT/m/lOOA 


13.41 mT/m/lOOA 


Gradient Strength at 500 A 


67.08 mT/m 


450 A 


60.37 mT/m 


330 A 


44.27 mT/m 


Number of Discrete Loops 


10 


Cu thick, between adjacent loops 


> 12.5mm 


Total Inductance (cable included) 


73nH 
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Total Resistance(Cable included) 


61mfl 


Rise Time at 2000V/500A 


19 jasec 


700V/450A 


49 (isec 


400V/330A 


64 |isec 


Linear Slew Rate at 2000V/500A 


3530 T/m/sec 


700V/450 A 


1235 T/m/sec 


400V/330A 


698 T/m/<iec 

V/-/ <J X / 111/ Owv* 


Duty Cycle (RMS Gradient) 




(500 A ref. point) at 4500 kW 


61% of the peak gradient(40.66mT/m) 


at 2000 kW 


40% of the peak gradient(27.1mT/m) 


at 1000 kW 


29% of the peak gradient(19.2mT/m) 



II: Table 6 

% 2) Shielded Biplanar Gradient Coils 



V| i)Transverse Y biplanar gradient coil : 

'4 For the shielded transverse biplanar gradient coils, the gap between the two primary 

coil's planes is chosen to be 2a=288mm, while the gap between the shielding planes is 

■| 6 2b=337mm. The constraint points for the minimization algorithm are shown in Table 7, The 

f if gradient field at the center of the structure is demanded to be 40 mT/m with 10% on axis 

j,!^ linearity at a distance of ± 130mm from the center of the gradient and 20% of f axis uniformity 

i!;!^ inside a 25 cm Diameter Spherical Volume (DSV). 



Constraints for a Shielded Biplanar Y gradient coil 


Xj in meters 


Yj in meters 


Zj in meters 


Bzscj in Tesla 


0.000 


0.001 


0.000 


0.0000400 


0.000 


0.130 


0.000 


0.0046800 


0.125 


0.001 


0.000 


0.0000320 


0.000 


0.001 


0.125 


0.0000320 



10 Table 7 

11 

1 2 The application of the Lagrange minimization technique to the set of constraints of the 

13 Table 7, generates a continuous current distribution for primary coil of the Transverse Shielded 

14 Y biplanar coil. With the assist of the stream function technique, the continuous current densities 
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1 for the primary coil can be approximated by a 15 discrete loops with a common current of 425.85 

2 Amps(figure 5). For a shielded Y biplanar gradient coil with parabolic returns, the current 

3 pattern for the primary coil is shown in Figure 6. while for the secondary coil its continuous 

4 current density can be approximated by 10 loops (figure 7b)with the same constant current as the 

5 primary coils. For a shielded Y biplanar gradient coil with parabolic returns, the current 

6 pattern for the secondary coil is shown in Figure 7b. Table 8 illustrates all the vital 

7 characteristics needed for the engineering and manufacturing phase of such a gradient coil. In 

8 addition figure 8 shows the on axis field behavior of the coil as it has been evaluated by applying 

9 the Biot-Savart law to the coil's discrete current pattern. 
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Shielded Y Biplanar Gradient Coil 




Property 




□ 


Gap primary/shielding 


288 mm/337mm 




Gradient Strength at mT/m/lOOA 


9.39 mT/m/lOOA 




Gradient Strength at 500 A 


46.95 mT/m 




450 A 


42.26 mT/m 




330 A 


31.00 mT/m 


-■y 


Number of Discrete Loops 


15/10 


mt 


Cu thick, between adjacent loops 


11.4mm/15.8mm 




Total Inductance (cable included) 


105^H 




Total Resistance(Cable included) 


95 mQ 




Rise Time at 2000V/500A 


27 ^tsec 




700V/450A 


72 i^sec 




400V/330A 


94 fisec 




Slew Rate at 2000V/500A 


1746 T/m/sec 




700V/450 A 


587 T/m/sec 




400V/330A 


330 T/m/sec 




Duty Cycle (RMS Gradient) 






(500 A ref. point) at 4500 kW 


47% of the peak gradient(22.27mT/m) 




at 2000 kW 


32% of the peak gradient(14.85mT/m) 
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at 1000 kW 


22% of the peak gradient(10.50mT/m) 


% Eddy Currents at 25cm DSV 


0.32% 



1 Tables 

2 

3 iv) Transverse X biplanar gradient coil: 

4 For the shielded transverse biplanar gradient coils, where the gap between the two 

5 primary planes is chosen to be 2a=284nnn, while the gap between the shielding planes is 

6 2b==341mm. The constraint points for the minimization algorithm are shown in Table 9. The 

7 gradient field at the center of the structure is demanded to be 40 mT/m with 10% on axis 

8 linearity at a distance of ± 130mm firom the center of the gradient and 20% off axis uniformity 

9 inside a 25 cm Diameter Spherical Volume (DSV). 

10 





Constraints for a Shielded Biplanar X gradient coil 




Xj in meters 


Yj in meters 


Zj in meters 


Bzscj in Tesla 


!«!? 

4 


0.001 


0.000 


0.000 


0.0000400 




0.130 


0.000 


0.000 


0.0046800 


:^ 


0.001 


0.125 


0.000 


0.0000320 


3 


0.001 


0.000 


0.125 


0.0000320 



li Table 9 



11 

li The application of the Lagrange minimization technique to the set of constraints of the 

14 Table , generates a continuous current distribution for primary coil of the Transverse Shielded X 

15 biplanar coil. With the assist of the stream function technique, the continuous current density for 

16 the primary coil can be approximated by a 16 discrete loops with a common current of 339.2 

17 Amps(figure 9). The continuous current density of the shielding coil can be also approximated 

18 by a set of 11 discrete loops carrying the same current as the primary coil(figure 10). Table 10 

19 illustrates all the vital characteristics needed for the engineering and manufacturing phase of 

20 such a gradient coil. 
21 

22 
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Shielded X Biplanar Gradient Coil 


Prooertv 




Gap primary/shielding 


284 mm/341nnn 


Gradient Strength at mT/m/lOOA 


11.79 mT/m/lOOA 


Gradient Strength at 500 A 


59.00 mT/m 


450 A 


53.10 mT/m 


330 A 


38.93 mT/m 


Number of Discrete Loops 


16/11 


Cu thick, between adjacent loops 


10.8mm/l 5.8mm 


Total Inductance (cable included) 


203(iH 


Total Resistance(Cable included) 


103 mQ 


Rise Time at 2000V/500A 


53 }j.sec 


700V/450A 


140 lasec 


400V/330A 


1 83 usee 


OICW JxclLC ctL Z-v/UvJ V / J\J\Jt\ 


1 1 %1 T/m/<ipr 

X 1 JL / lil/ oCL- 


70nV/4^0 A 

/ UU V / *-tJ\J t\ 


'^jLjL X / XIX/ dvv 


400V/330A 


213 T/m/sec 


Duty Cycle (RMS Gradient) 




(500 A ref. point) at 4500 kW 


45% of the peak gradient(26.55mT/m) 


at 2000 kW 


30% of the peak gradient(17.59mT/m) 


atlOOOkW 


21% of the peak gradient(12.44mT/m) 


% Eddy Currents at 25cm DSV 


0.37% 



Table 10 



v) Axial Z biplanar gradient coil: 

For the shielded Axial biplanar gradient coils, where the gap between the two primary 
planes is chosen to be 2a=292mm. The distance between the two shielding planes is 2b=345mm. 
The constraint points for the minimization algorithm are shown in Table 1 1 . The gradient field at 
the center of the structure is demanded to be 40 mT/m with 10% on axis hnearity at a distance of 



28 



1 +1 SOinm from the center of the gradient and 7.5% of f axis uniformity inside a 25 cm Diameter 

2 Spherical Volume (DSV). 
3 



Constraints for a Shielded Biplanar Z gradient coil 


Xj in meters 


Yj in meters 


Zj in meters 


Bzscj in Tesla 


0.000 


0.000 


0.001 


0.0000400 


0.000 


0.000 


0.130 


0.0046800 


0.125 


0.000 


0.001 


0.0000394 


0.000 


0.125 


0.001 


0.0000394 



4 Table 11 

5 

The application of the Lagrange minimization technique to the set of constraints of the 
Table 1 1, generates a continuous current distribution for primary coil of the Axial Shielded Z 
fM biplanar coil With the assist of the stream function technique, the continuous current densities 
I j for the primary coil can be approximated by 1 6 discrete loops with a conmion current of 420.9 
•lb Amps(figure 11). the continuous current density for the shielding coil can be approximated by 12 
:il 1 discrete loops (figure 12) Table 12 illustrates all the vital characteristics needed for the 
JJZ engineering and manufacturing phase of such a gradient coil. 
'% 

m 
i'i4 



Shielded Z Biplanar Gradient Coil 


Property 




Gap primary/shielding 


292 mm/345mm 


Gradient Strength at mT/m/lOOA 


9.50mT/m/100A 


Gradioit Strength at 500 A 


47.51 mT/m 


450 A 


42.76 mT/m 


330 A 


3 1.36 mT/m 


Number of Discrete Positive Ix)ops 


16/12 


Cu thick, between adjacent loops 


ILOmm/l 5.8mm 


Total Inductance (cable included) 


108 nH 
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Total Resistance(Cable included) 


107 mQ 


Rise Time at 2000V/500A 


28 jisec 


700V/450A 


75 usee 


400V/330A 


98 |j.sec 


Slew Rate at 2000V/500A 


1712 T/m/sec 


700V/450A 


573 T/m/sec 


400V/330A 


320 T/m/sec 


Duty Cycle (RMS Gradient) 




(500 A ref. point) at 4500 kW 


44% of the peak gradient(20.68mT/m) 


at 2000 kW 


29% of the peak gradient(13.79mT/m) 


at 1000 kW 


21% of the peak gradient(9.75mT/m) 


% Eddy Currents at 25cm DSV 


0.33% 



Table 12 



3) Non-Shielded Uniplanar Gradient Coils 

i) Transverse X, Y uniplanar gradient coil 

For the non-shielded transverse uniplanar gradient coils, where the gap between the plane 
from the geometric center of the structure of the gradient coil is chosen to be -a —6mm, the 
constraint points for the minimization algorithm are shown in Table 13. The gradient field at the 
center of the structure is demanded to be 40 mT/m with 10% on axis linearity at a distance of ± 
100mm from the center of the gradient and 20% of f axis uniformity inside a 25 cm Diameter 
Spherical Volume (DSV). 



Constraints for a Non-Shielded Uniplanar transverse(X,Y) gradient coils 


Xj in meters 


Yj in meters 


Zj in meters 


Bzscj in Tesla 


0.000 


0.001 


0.100 


0.0000400 


0.000 


0.060 


0.100 


0.0024000 


0.000 


0.200 


0.100 


0.0072800 


0.180 


0.001 


0.100 


0.0000360 
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0.000 


0.001 


0.150 


0.0000410 


0.000 


0.001 


0.200 


0.0000350 



1 Table 13 

2 

3 The application of the Lagrange minimization technique to the set of constraints of the 

4 Table 13, generates a continuous current distribution for primary coil of the Transverse Non- 

5 Shielded Y uniplanar coil. With the assist of the stream function technique, the continuous 

6 current densities for the coil can be approximated by a 10 discrete loops with a common current 

7 of 41 1.2 Amps (figure 13). Table 14 illustrates all the vital characteristics needed for the 

8 engineering and manufacturing phase of such a gradient coil. In addition figure 14 shows the on 

9 axis field behavior of the coil as it has been evaluated by applying the Biot-Savart law to the 
1 0 coil's discrete current pattern. 

o 





Non-Shielded Transverse (X,Y) Uniplanar Gradient Coil 


""''■4 


Property 












Gradient Strength at mT/m/lOOA 


9.73 mT/m/lOOA 




Gradient Strength at 500 A 


48.64 mT/m 


iU 


450 A 


43.77 mT/m 


i 


330 A 


32.11 mT/m 


!■• ; 


Number of Discrete Loops 


10 




Cu thick, between adjacent loops 


> 10.2mm 




Total Inductance (cable included) 


295m,H 




Total Resistance(Cable included) 


73 m^^ 




Rise Time at 2000V/500A 


75 (jsec 




700V/450A 


199 usee 




400V/330A 


259 usee 




Slew Rate at 2000V/500A 


648 T/m/sec 




700V/450A 


220 T/m/sec 




400V/330A 


124 T/m/sec 
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Duty Cycle (RMS Gradient) 




(500 A ref. point) at 4500 kW 


53.4% of the peak gradient(26 mT/m) 


at 2000 kW 


36% of the peak gradient(17.33mT/m) 


at 1000 kW 


25% of the peak gradient(12.26mT/m) 


Table 14 



2 

3 ii) Axial Z Uniplanar gradient coil: 

4 For the non-shielded Axial Uniplanar gradient coil, the gap between the geometric center 

5 of the structure of the gradient coil and the planar surface is chosen to he -a'=-6mm. The 

6 constraint points for the minimization algorithm are shown in Table 15. The gradient field at the 

7 center of the structure is demanded to be 40 mT/m with 20% on axis linearity at a distance of + 

8 200mm Jfrom the center of the gradient with no-rollover point and 20% of f axis uniformity 
% inside a 20 cm Diameter Spherical Volume (DSV). 

M 





Constraints for a Non-Shielded Uniplanar Z gradient coil 




Xj in meters 


Yj in meters 


Zj in meters 


Bzscj in Tesla 




0.000 


0.000 


0.031 


0.0000400 


ill 


0.000 


0.000 


0.231 


0.0062800 




0.180 


0.000 


0.031 


0.00003600 




0.000 


0.200 


0.031 


0.00003200 



12 

13 The application of the Lagrange minimization technique to the set of constraints of the 

14 Table 15 generates a continuous current distribution for primary coil of the Axial Non-Shielded 

1 5 Z uniplanar coil. With the assist of the stream function technique, the continuous current 

16 densities for the two coils can be approximated by a 10 discrete loops with a common current of 

17 316.1 Amps(figure 15). Table 16 illustrates all the vital characteristics needed for the 

1 8 engineering ^d m^ufacturing phase of such a gradient coil. 
19 
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Non-Shielded Z Uniplanar Gradient Coil 


Property 








Gradient Strength at mT/m/lOOA 


12.65 mT/m/lOOA 


Gradient Strength at 500 A 
450 A 
330 A 


63.25 ml/m 
56.93 mT/m 
41.75 mT/m 


Number of Discrete Loops 


13 


Cu thick, between adjacent loops 


> 12.5mm 


Total Inductance (cable included) 


241 |iH 


Total Resistance(Cable included) 


57 mD 


Rise Time at 2000V/500A 
700V/450A 
400V/330A 


62 fjsec 
160 usee 
209 fjsec 


Linear Slew Rate at 2000V/500A 
700V/450 A 
400V/330A 


1034 T/m/sec 
354 T/m/sec 
200 T/m/sec 


Duty Cycle (RMS Gradient) 
(500 A ref. point) at 4500 kW 
at 2000 kW 
at 1000 kW 


63% of the peak gradient(40.0mT/m) 
42% of the peak gradient(26.7mT/m) 
30% of the peak gradient(18.9mT/m) 



Table 16 

4) Shielded Uniplanar Gradient Coils 



i)Transverse (X,Y) uniplanar gradient coil with parabolic cut-offs: 

For the shielded transverse biplanar gradient coils, the gap between the geometric center 
of the entire gradient structure and the planar surface of the primary coil is chosen to be -a'=- 
6nmi. The gap between the geometric center of the structure and the planar surface of the 
shielding coil is chosen to be -b'=-70mm. The constraint points for the minimization algorithm 
are shown in Table 17. The gradient field at the center of the structure is demanded to be 40 
mT/m with 20% on axis linearity at a distance of 
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1 + 200mm from the center of the gradient. In addition, the off-axis uniformity of the gradient field 

2 is restricted to within 20% from the field's ideal value inside a 25 cm Diameter Spherical 

3 Volume (DSV). 



Constraints for a Shielded Uniplanar transverse(X,Y) gradient coils 


Xj in meters 


Yj in meters 


Zj in meters 


Bzscj in Tesla 


0.000 


0.001 


0.075 


0.0000400 


0.000 


0.100 


0.075 


0.0040000 


0.000 


0.300 


0.075 


0.0140400 


0.180 


0.001 


0.075 


0.0000380 


0.000 


0.001 


0.13 


0.0000390 


0.000 


0.001 


0.180 


0.0000320 



4 Table 17 

9 

'f0 The application of the Lagrange minimization technique to the set of constraints of the 

J Table 17 generates a continuous current distribution for primary coil of the Transverse Shielded 

biB Y biplanar coil. With the assist of the stream fimction technique, the continuous current densities 

I ^ for the primary coil can be approximated by a 1 6 discrete loops with a common current of 324, 1 7 

::H) Amps(figure 16a), while for the secondary coil its continuous current density can be 

fP approximated by 8 loops (figure 16b)with the same constant current as the primary coils. Table 
18 illustrates all the vital characteristics needed for the engineering and manufacturing phase of 

i such a gradient coil. In addition figure 16c shows the on axis field behavior of the coil as it has 

14 been evaluated by applying the Biot-Savart law to the coil's discrete current pattern. 

15 



Shielded Transverse (X,Y) Uniplanar Gradient Coil 


Property 








Gradient Strength at mT/m/lOOA 


12.33 mT/m/lOOA 


Gradient Strength at 500 A 
450 A 
330 A 


61.69 mT/m 
55.52 mT/m 

40.70 mT/m 
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Number of Discrete Loops 


16/8 




Cu thick, between adjacent loops 


12.4mm/15.8mm 




Total Inductance (cable included) 


302pH 




Total Resistance(Cable included) 


122 mQ 




Rise Time at 2000V/500A 


78 usee 




700V/450A 


210 usee 




400V/330A 


277 usee 




Slew Rate at 2000V/500A 


792 T/m/sec 




700V/450 A 


264 T/m/sec 




400V/330A 


147 T/m/sec 




Duty Cycle (RMS Gradient) 






(500 A ref. point) at 4500 kW 


40% of the peak gradient(24.95mT/m) 




at 2000 kW 


27% of the peak gradient(16.64mT/m) 


:' \9i 


at 1000 kW 


19% of the peak gradient(11.76mT/m) 




% Eddy Currents at 25cm DSV 


0.39% 


fi Table 18 

a 




1 !.U 

jll iii) Axial Z biplanar gradient coil: 



i:| For the shielded Axial Uniplanar gradient coils, the gap between the geometric center of 

^§ the entire gradient structure and the planar surface of the primary coil is chosen to be -a'— 6mm. 

7 The gap between the geometric center of the structure and the planar surface of the shielding coil 

8 is chosen to be -b --70nmi. The constraint points for the minimization algorithm are shown in 

9 Table 19. The gradient field at the center of the structure is demanded to be 40 mT/m with 20% 

10 on axis linearity at a distance of ± 1 30mm from the center of the gradient and 7.5% of f axis 

1 1 uniformity inside a 25 cm Diameter Spherical Volume (DSV). 
12 



Constraints for a Shielded Uniplanar Z gradient coil 


Xj in meters 


Yj in meters 


Zj in meters 


Bzscj in Tesla 


0.000 


0.000 


0.031 


0.0000400 
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0.000 


0.000 


0.231 


0.0062800 


0.180 


0.000 


0.031 


0.00003600 


0.000 


0.200 


0.031 


0.00003200 



1 Table 19 

2 

3 The application of the Lagrange minimization technique to the set of constraints of the 

4 Table 19, generates a continuous current distribution for primary coil of the Axial Shielded Z 

5 biplanar coil. With the assist of the stream function technique, the continuous current densities 

6 for the primary coil can be approximated by 21 discrete loops with a common current of 385.71 

7 Amps(figure 1 7a). The continuous cuirent density for the shielding coil can be approximated by 

8 13 discrete loops (figure 17b) Table 20 illustrates all the vital characteristics needed for the 
, 9 engineering and manufacturing phase of such a gradient coil. 

m 



Shielded Z Uniplanar Gradient Coil 


Property 








Gradient Strength at mT/m/lOOA 


10.37 mT/m/lOOA 


Gradient Strength at 500 A 
450 A 
330 A 


51.85 mT/m 
46.67 mT/m 
34.22 mT/m 


Number of Discrete Positive Loops 


21/13 


Cu thick, between adjacent loops 


9.4nmi/15.8nun 


Total Inductance (cable included) 


306 nH 


Total Resistance(Cable included) 


125 mQ 


Rise Time at 2000V/500A 
700V/450A 
400V/330A 


80 |Lisec 
213 (Asec 
281 ^sec 


Slew Rate at 2000V/500A 
700V/450 A 


656 T/m/sec 
218 T/m/sec 
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400V/330A 


122 T/m/sec 


Duty Cycle (RMS Gradieait) 
(500 A ref. point) at 4500 kW 
at 2000 kW 

at 1000 kW 


40% of the peak gradient(20.97mT/m) 
27% of the peak gradient(13.98mT/m) 
19% of the peak gradient(9.89mT/m) 


% Eddy Currents at 25cm DSV 


0.19% 



1 Table 12 

2 

3 Mutual Inductance Calculations 
4 

5 Bivlanar Gradients 
6 

% The mutual energy between two biplanar gradients of the same axis was evaluated based 

■if on the equations (64) and (65). In particular for this example the mutual energy between two 

i2 biplanar transverse Y gradients was evahiated. For the first biplanar set, the gap between its 

lOi primary planes was set at 2a=288mm, while the gap between its shielding planes was set at 

W 2b=341mm. For the second biplanar set, the gap between the primary planes was set at 

ll 2a=282mm, while the gap between the two secondary planes is set at 2 b =292mm. The mutual 

1| energy of these two sets as a function of shifting one set with respect to the other along the z 

i i direction. The optimum (zero) mutual energy of these two sets occurs at the point zo=216mm 

1| from the geometric center of the fir gradient coil set. 
16 

17 Uniplanar Gradients 
18 

19 The mutual energy between two uniplanar gradients of the same axis was evaluated based 

20 on the equations (21). In particular for this example the mutual energy between two uniplanar 

21 axial Z gradients was evaluated. For the first uniplanar set, the gap on the primary plane was set 

22 at -a —6mm, while the gap of the shielding plane was set at -b --70mm. For the second biplanar 

23 set, the gap for the primary plane was set at -a —3mm, while the gap between the two secondary 

24 planes is set at - & — 75nmi. The optimum (zero) mutual inductance for these two set occurs at 

25 the point zo=245mm from the geometric center of the first gradient coil set. 
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1 

2 The present invention describes a method of designing and manufacturing an Open 

3 Architecture gradient coil set which is suitable for magnetic resonance systems with vertically 

4 directed or horizontally directed fields. The concept of designing of an open architecture gradient 

5 coil system offer a significant advantage over existing geometries since, the present gradient 

6 design provides superior performance in terms of gradient strength and field linearity and 

7 uniformity while the structure of the coil is an open architecture focusing of reducing patient 

8 claustrophobia. The invention can be implement in gradient coil designs for open magnet 

9 geometries with either vertically and or horizontally directed fields. In particular, for magnets 

10 with vertically directed fields, besides the open architecture design of the proposed gradient 

1 1 system a significant gradient increase in the strength of the local gradient field up to 55mT/m 

12 with slew rates exceeding 1000 T/m/sec in an imaging volume covering the size of the human 
tl head (25 cm DSV). Another significant advantage of the design is the abiUty to control and 

W minimize the dB/dt and eddy current levels for shielded design. In addition, the open architecture 

i S design can be paired in a phased array configuration with altemative shielded or not open 

1 1 architecture gradient coil set in such a configuration that covers the human body head to toe. The 

significant advantage of such configuration is the significant reduction of dB/dt effects and the 

|§ utilization of high strength non-claustrophobic gradient coils for high speed imaging of the 

II human body with excellent gradient field uniformity and linearity characteristics in the local 

in predetermined gradient volume. For open magnets with horizontally directed fields (double- 

^1 donut designs), the invention provides a competitive advantage over existing gradient designs 

22 since it can deliver gradient fields with strengths two to three fold grater than akeady existing 

23 designs and the ability to cover the entire human body in a phased array configuration with 

24 minimized dB/dt effects and induced eddy currents on the magnet's shields or metallic surfaces 

25 during the locahzation of the gradient field. The open architecture can be modified with 

26 paraboKc cut-offs in order to incorporate the human shoulders for the portion of the open 

27 architecture gradient system which is suitable for imaging the human head. Furthermore, the 

28 design can be incorporated as an insertable configuration and can be fully integrated with any 

29 volumetric or surface phased array or singular coil configuration for imaging portions or the 

30 entire human body. In another aspect of the invention the design can be secured into the structure 

31 of its corresponding design. Although the design incorporates a limited curvature for the 
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1 perspective gradient coil configurations, the curvature can be easily modified for those skilled to 

2 tiie art to any unspecified shape suitable for the any targeted application. In addition, for a self- 

3 shielded design the side planes of the secondary (shielded) coil can be mechanically shifted and 

4 can be adjusted with respect the primary planes of the coil in an array configuration with zero 

5 mutual inductance in such a way that the imaging volume of the coil is extended. 

6 The present invention provides a gradient coil set that minimizes patient anxiety while 

7 delivering high peak, uniformity and slew rate gradient coils for Large imaging volumes. The 

8 combination of both biplanar and uniplanar gradients achieve a gradient magnetic field with 

9 superior field qualities. The coil set may have the side planes with parabolic cut offs in order to 

10 assist in patient comfort for applications that are focused on imaging the head, the neck, the 

1 1 spine, the heart and the upper torso area. Due to the phased array design low dB/dt levels can be 

12 achieved. In the gradient coil set and especially for the shielded designs, the outer shielding 

il planes can be moved to form a phased array configuration with their perspective primary planes 

11 corresponding to the same gradient coil axis. The movement can be done via mechanical means 

i:i or any other means. An unlimited spectrum of FMRI and Literventional applications can be 

I | performed with the he "Open Face Architecture" gradient coil, because of its openness and the 

II high peak and slew rate values for the gradient field. 



Referring to FIG. 1 , a gradient coil set 10 includes a uniplanar Z-gradient coil 12 and a 

il biplanar X-gradient coil 14 and a biplanar Y-gradient coil 16. Each of coils may include a 

iS) respective shield coil 1 8, 20, 22. The coil set 10 has an open Z-axis face opposite the coil 12. 

IJ The biplanar coils may be provided with shoulder reUefs 24. The structure of the biplanar coils 

22 14, 16 may be conjoined with the structure of the uniplanar coil 12, for example, as shown. 

23 Referring to FIG. 1 8, a gradient coil set 30 includes self-shielded gradient coils 32, 34, 36 

24 and a RF head coil 38 

25 Referring to FIG. 1 9, a gradient coil set 40 includes un-shielded gradient coils 42, 44, 46 

26 and a RF head coil 48. 

27 Referring to FIG. 20, a gradient coil set 50 includes self-shielded gradient coils 52, 54, 

28 56, a RF head coil 58 and a CTL spine coil 59 integrated into the structure of the gradient coil set 

29 50. 

30 Referring to FIG. 21, a gradient coil set 60 includes un-shielded gradient coils 62, 64, 66, 

31 a RF head coil 68 and a CTL spine coil 69 integrated into the structure of the gradient coil set 60. 
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Referring to FIG. 22, a gradient coil set 70 includes Z-gradient uniplanar coils 72, 74 and 
X-Y biplanar gradient coils 76, 78 in a phased array configuration. A head RF coil 79 may be 
included. 

Referring to FIG. 23, an open Z-axis face coil set 90 includes a knee and upper thigh RF 
coil 92 integrated therewith. 

Referring to FIG. 24, an open Z-axis face coil set 100 includes a knee and foot RF coil 

102 integrated therewith. 

It should be evident that this disclosure is by way of example and that various changes 
may be made by adding, modifying or eUminating details without departing from the fair scope 
of the teaching contained in this disclosure. The invention is therefore not limited to particular 
details of this disclosure except to the extent that the following claims are necessarily so limited. 
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